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ABSTRACT 


Empirical evidence is presented that solar wind thermal electrons obey 

a polytrope law with polytrope index y = 1-175 ± 0.03 (3o). The Voyager 2 

and Mariner 10 data used span the radial range from 0.^15 AU to ^.?6 AU and 

have a large dynamic range in density (4 decades) , and over one decade in 

temperature, which is crucial for an unambigious determination of y • We 

find that there is large xt' ± 50% variation in the entropy like stream 

tube constant 2 relative to its best fit value For the Voyager 2 data 

set, which has 5831 hourly averages, this variation comprises a nearly 

unbiased statistical sample of log S relative to log Z^, thereby reducing 

systematic errors . A positive correlation on average between E and bulk 

speed is found, which does introduce systematic errors (v/' 5%) in the 

determination of y resulting in larger uncertainties than implied by a 

simple regression of log T vs log n. Our result supports the thei;)retical 

€ 

predictions by Scudder and Olbert (1979b) that solar wind thermal electrons 
in the asymptotic solar wind, should obey a polytrope law with polytrope 
index y = 1.l5 ± 0.06 (3o) . It does not support the results reported by 
Feldman ei (1978) for high speed streams where they find 

y = 1,45 ± 0,15 (3°). We attribute this difference to the biased sample of 
stream tube constants 2 encountered in that limited study and/or the 
smaller dynamic range of density and temperature available to the 1 AU 
observer. Our results support; 1) the widespread impressions in the 
literature that solar wind electrons behave more like an isothermal than 
adiabatic gas; 2) the arguments put forward by Scudder and Olbert (1979ajb) 
that Coulomb collisions are the dominant stochastic process shaping 
observed electron distribution functions in the solar wind, and 3) the 
assignment if the interplanetary potential as equal to approximately seven 
times the temperature of the thermal electrons , 
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I. INTRODUCTION 

Attempts to model the two-fluid solar wind expansions have either used 
the full electron and ion energy equation with Spitzer's formulation of 
thermal conduction (Hartle and Sturrock (1966)) or have simplified the 
energy equations by assuming two different polytropic laws for ions and 
electrons (Goldstein and Jokipii (197?) and V. J. Pizzo, private 
communication). Both of these approaches have fundamental difficulties. 
Concurrent with the existence of a Spitzer thermal conductivity is the 
necessary condition that the Coulomb free path over the scale length be 
much less than unity and that the system be slightly removed from 
(homogeneous) thai raodynamic equilibrium with local Maxwellian velocity 
distribution functions; however, the observed distribution functions do not 
have this property. The second approach represents a considerable 
mathematical simplification of the fluid equations; however, they must be 
regarded as phenomenological. In addition, since any heat law can be 
imitated by a spatially varying polytropic index, Y, the use of a single 
constant y represents an additional assumption of the specific form of the 
heat law, namely that the heat transport is implemented as an enthalpy flux 
(cf . Parker, 1963) • 

In this paper we present empirical evidence that between 0.4 and 5 AU 
the thermal portion (but not all ) of the solar wind electron population 
obeys a polytrope relation. We also show that this functional relationship 
is a member of a broader class of possible laws required of a steady state, 
fully ionized plasma whose proper frame electric field is dominated by the 
polarization electric field. The empirically determined, thermodynamically 
interesting value of Y - 1.175 is virtually that predicted (y = 1.16) hy 
the theoretical considerations of Scudder and Olbert (1979b). We thus 
provide strong, direct, empirical evidence for the nearly isothermal 
behavior of solar wind electrons as has been indirectly argued in the 
literature for some time (Burlaga^^. (1971), Hundhausen and Montgomery 
(1971), Hundhausen (1973))- 

O' 

This polytrope relation for the thermal electrons refers to the bulk 
95^) of the number density and 7556 of the total eldctron pressure. By 
contrast the heat flow vector requires an accurate description of all the 


electrons including that peculiar minority of particles which make the 
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electron velocity distributions function skewed in the proper frame. This 
group of particles is usually in the trans to extrathermal regime C > kT: 
(observationallv ktontgomerv et al . (1968)^ Scudder (1970), Montgomery 
(1972); theoretically Scudder and Olbert , (1979a), hereafer referred to as 
SO-I). For the structureless solar wind in the regime where conduction is 
not energetically important the electron thermodynamics may in first 
approximation be described by this polytrope relation. However, the 
detailed description of the initial coronal acceleration or dynamical 
interactions in structured flow at any distance requires a further 
elaboration of the actual heat law. In particular the inference that 
advection is the principle mode of internal energy transport for electrons 
in stream interaction regimes beyond 20-30 is probably inappropriate. 
This may be seen since use of the polytrope law stream line by stream line, 
neglects the possibility that heat can be locally deposited along the 
stream tube in the form of heat without the correlative amount of density 
enhancement . 


2. THEORETICAL BACKGROUND 

Because electrons have a small inertia relative to ions , the proper 
frame electric field present in the steady state inhomogeneous solar wind 
expansion is inextricably connected to the thermodynamic forces upon these 
particles. The generalized Ohm's law is the governing relation between 
this electric field, E* , in the proper frame and the gradient of 
macroscopic variables which are the local thermodynamic forces. The 
relation is determined by the difference of the electron and ion momentum 
equations and in general is given by equation 1 
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where, E and E* are the electric field in the inertial and proper frame, 
respectively, 5 the magnetic field, u the proper frame velocity, the 
electron pressure tensor, n the electron density, j* the proper frame 

® -f 

current density, n the free charge density, g the gravitational force per 
unit mass, e the unit electric charge, m^ the 
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eleobron mass, and = + u ’ V is the total time derivative. 

dt 3t 

The terms on the RHS of Eq. 1 have been written in decreasing order of 
usual importance in the steady state solar wind as argued in Rossi and 
Olbert ( 1970 ), where most notably the first three terms are, respectively, 
the polarization. Hall, and collisional (or Olimic heating) terras. The 
remaining set of terms are usually negligible in astrophysical 
circumstances . 

The collisional group may be dropped since it is inversely proportional 
bo the magnetic Reynolds number which is large in the solar wind. The size 
of the Hall terra has been estimated assuming Alfvenic fluctuations and 
shown to be small 10^) relative bo the polarization berm except, 
possibly in the lower corona. 

In this section we will explore the implications of E* being 
predominately determined by the polarization term. 

S* = -V • ?„/en„ (2) 

Using an argument made by S. Olbert (private communication, 1978) the curl 
free nature of E* in steady state with P = Pg I and Eq . 2 implies that 

(-VP ) (-VP ) 1 

VxE» = 0«Vx ^ Vx ^ = + (VPg) X (vn^) (3) 

en en en ^ 

a e e 

Therefore, P^ = f(Hg) is the most general relation which satisfies this 
requirement. A particular example of such a function f(n^) is the 
polytrope law 

Pg = S nj (4) 

where 2 is a stream tube constant and Y the polytrope index. However, 
without further justification there appears no rationale other than 
simplicity that a polytrope relation such as (4) should be the particular 
function f observed to relate thermal electron temperature and ambient 
density which satisfies (3) . 

The properties of the velocity dependence of Coulomb scattering were 
examined in Scudder and Olbert (1979b), hereafter referred to as SO-II 
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where it was shown both theoretically and empirically that a change in the 
solar wind electron velocity distribution function occurs in the vicinity 
of a transition energy = 7 kT„ , where T is the "temperature" of the 
thermal (core) electrons; the constant J only depended on fundamental 
constants (m„, m„, e, and that T_) . These authors explored the 

6 p Q P 

implications that this dime)r-’'ionless transition energy g»/kT^ could be 
independent of space and timv* wnile was simultaneously the 
interplanetary potential suggested by exosphere approaches to the solar 
wind electron theory: (Jockers (1970); Schultz and Eviatar (1972); Perkins 

(1973)) and observations: (Feldman (1975), Rosenbauer ^ 

( 1976 )). It was concluded in SOU that the necessary and sufficient 
condition that this be possible is that the thermal electrons do obey a 
polytrope law, Eq. 4, with an expected Y = 1.16 ± O.O 6 (3cr) . We now 
proceed to test their conclusions. 

3 . SOLAR WIND ELECTRON OBSERVATIONS VOYAGER 2 - MARINER 10 

The simplest manner of empirically determining y would be a linear 
least squares fit of the two observables log T^ versus log n^ with best 
slope, and its uncertainties being identified with y by 


Y 


“fit * ’ 


(5) 


Because the variation of T^ with radius previously reported by Gringauz and 
Verigin (1975) Ogilvie and Scudder (1978) and Sittler jet ( 1979) is so 
shallow, obtaining a significant dynamic range in temperature (greater than 
one order of magnitude) requires a protracted period of data collection for 
a deep space mission such as Voyager. The Voyager- results to be discussed 
have been derived from data collected by a cylindrical Faraday cup between 
1.38 and 4.76 AU over a period of 1 I /3 years, or about I 8 solar rotations. 
These data span one order of magnitude variation in thermal temperature and 
three orders of magnitude in density. The Mariner 10 data were collected 
by a spherical section electrostatic analyzer between 0.9 and 0.45 AU over 
a period of 2 months . These data extend the composite data coverage to 
nearly four orders of magnitude in density and over an order of magnitude 
in thermal electron temperature. (Further details of Voyager and 
Mariner 10 processing and experimental assumptions are described in 
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(Sibtler and Scudder (1979) and Ogilvie and Scudder (1978); the Voyager 
instrument is described in detail in Bridge ^ , 1977))- 

Because neither of the variablee in .the suggested relation (4). are , 
detemirjed with absolute precision, the linear least squares approach for 
determining y is inadequate. The more general .approach of "hypothesis 
testing" (Eadie jgt, al . , (1971)) is required whereby the perpendicular 
dispersioq of the data from the optimal hypothesis is minimized. This 
procedure corresponds in this case to determining the polytrope constant z 
(intercept) and index y>=. -( slope + 1) in a way that minimizes the 
perpendicular dispersion of our observations relative to this hypothesized 
law. In log-log space the polytrppe hypothesis along eaph stream tube 
labeled "s". is a straight line 

log T(s.) = log n (s) + lo^ £ (s), . . (6a) 


= “fit ^ 


(6b) 


with slope and intercept I. Minimizing the perpendicular dispersion 

about this .hypothesis is particularly simple. 


Voyager 2 Data 

We have performed this minimizing procedure on the Voyager data and the 
results are shovjn in Figure 1 . In the body of the figure is plotted the 
optimal polytrope hypothesis with empirical polytrope index y = Ofj^t ■ " 
1.185 i 0.018 (3<J) and with optimal average stream tube constant 
< log Z > = I given by M.74 ± 0.009 (3d) where the units of Z are 
determined . by T (°K) and n (cm ),. Each cross represents the mean 
determined from the Voyager 2 hourly averaged data points which are found 
in bins perpendicular to and of equal width along the best, fit line. The 
583.1 Voyager, hourly data points were used to determine the optimal 
polytrope relation. As a group the binned data points and error bars (in 
many . cases .error bars are smaller than the size of the crosses) closely 
resemble the .polytrope relation selected by the above described 
minimization process . The error flags are parallel to the sides .pf the 
bins which the optimal hypothesis slope defines, being perpendicular to the 
selected best hypothesis line. (These lines do not appear perpendicular 
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since decades of density and temperature do not have the same physical 
dimension in the figure). It is, however, clear that attempts 'to determine 
the best fit line for subsets of the data with smaller than the current 
dynamic range could determine slopes that are grossly inconsistent with the 
trend of the composite set (of. Feldman ^ , 197&). 


Mariner 10 Data 

The solid dots in Figure 1 show the placement of all similarly binned 
hourly averaged Mariner 10 electron data in this log T, I'-*' n space. The 
Mariner 10 data taken in 1973 had a slightly higher average stream tube 
constant = 4.83 than for the average Voyager data = ^*7^ ± 

009 3<^). The Mariner 10 data have been placed on this figure so that as a 
set they have the same average stream tube constant as obtained by the 

XT' 

Voyager fit. This has been done by a slight 1/f) displacements of the 
direct observations perpendicular to the average Voyager polytrope fit 
relation. The Mariner 10 data points clearly follow the trend suggested by 
the Voyager data set and give further dynamic range to the empirical 
relation . 

In inset 1A of Figure 1 we address the uncertainty with which the 
polytrope constant < log 2 > and index y = + 1 have been determined. 

The panel is comprised of a gray scale map of the surface in the 
vicinity of the optimal < log 2 > and selected by the non-linear least 

squares procedure outlined above. Motion parallel to the horizontal axis 
corresponds to changing the slope of the hypothesis while keeping the 
intercept fixed. Motion along the ordinate corresponds to changing the 
intercept while keeping the slope fixed. The gray coded intensity 
corresponds to the magnitude of x^, the normalized mean square dispersion 
of the data about a polytrope hypothesis, with polytrope constant and index 
determined by its location on the X^ map. The variation of X^ in 
increments of intensity from 0 to 15 is given by 
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< X 
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where = 7.0? and = 8.08. Using the procedure discussed by 

Lampton ^ al . (1976), the 99^ confidence level occurs just inside the ± 3 '^ 
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range of the plot (where o has been defined by the linear least squares fit 
to the binned data. It is particularly important to note that there are no 
topological "channels" in the surface that lead away from the optimally 
selected and I which are found to be at the bottom of a nearly 

circular well in the dispersion surface. This is not a foregone conclusion 
because of the non-linear character of the minimization process. 

It is appropriate to document the distribution of the Voyager 

observations in the T, n space as available from the spacecraft time 

series . This has been graphically done in inset B of Figure 1 . The 

central square panel of inset IB depicts the intensity keyed probability 

(D. ./cr.) of departure of individual hourly averaged observations from 

the optimal polytrope relation shown in the main figure, a. is the 

variance of all the displacements D. . in the jth bin. The horizontal axis 

(or bin index j) is measured in units along the best fit line of the main 

figure. For every cross in the main figure there corresponds 3 vertical 

columns (bins) in the central and lower panel of inset IB. All 

observations were found to be confined to ± 3a; of the optimal position. 

The intensity keyed probability distribution P. in each column (bin) 

u 

represents the observed normalized frequency of departure within each bin 
of the stream tube constant 2(t) sampled by Voyager relative to the average 
< Z(t) > determined from the composite fit. 

The upper horizontal panel gives the 16 linearly graded levels of 
Intensity for the probability of occurrence. The lower horizontal panel 
gives the percent of the total Voyager data set within each bin, where the 
blank keys represent less than 29 observations within each bin while 
typically there are 170 observations within each bin. 

The vertical steps to the right of the central panel depicts in normal 
histogram form the composite distribution of stream tubes constants 
sampled . This composite is defined as 


1 


P (D /a) = 


N, 


r 


"j 


( 6 ) 


j=1 


where N™ is the total number of observations, n. the total number 

X J 

of points in each bin , and Ng the number of bins . 
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Inset Figure- IB documents the followihg; (!) {in the central panel) 
that the most aberrant points from the poly'tfope law are from data obtained 
over short enough time to have a biased set of stream tube constants for 
the regime encountered (the first three points of main Figure 1 come from a 
sparse high velocity stream at 4,76 AU at the very end of our time series); 
these data are thus strongly coherent and not statistically representative 
of stream tube constants for this density regime; (2) (lower strip) these 
biased data points are not dominant in determining our optimal polytrope; 

(3) adjacent subsets of bins over a small dynamic range may or may not have 
gaussian Pj (D^j/Oj) , i.e., impartial samples of stream tube constants; and 

(4) (right strip) that the entire data set as a whole represents a nearly 

statistically unbiased sample of stream tube constant regimes (i.e., note 
gaussian shape of composite P(Dj^/cr) where P(1) 1/e). To illustrate this 

result we have compiled in Figure 2 a plot of the percent of occurrence of 

““fit 

log 2 as determined by hourly averaged log (Tj^/(n^ ) with mean value 


< log 2 > = 4.74. The near gaussian shape with Icj width 50?S of < 2 > 
shown here, means the data set is comprised of a nearly statistically 
random distribution of log 2. However, as the equal variance reference 
gaussian curvvS of the figure illustrates, the sample does have a slight 
shewness in the sample frequencies of log 2 . 

The fourth conclusion above and Figure 2 suggest that the composite 
Voyager data represent a nearly statist^'^al sample of stream tube constants 
with respect to the "typical" best fit line which gives the average stream 
tube constant = 10 . Before the mathematically defined best fit slope, 
“fit’ interpretable as the thermodynamically interesting, physical, 
polytrope index along a given stream tube, we must quantify the degree of 
correlation present in the data between stream tube constants sampled and 
the density. 


4, FITS AND PHYSICS 

We must now focus on the interpretation of the fits presented in 
relation to our assumed model. The fits pertain to data collected across 
many different stream tubes in the solar wind, whereas the hypothesis being 
tested pertains to the physics of thermal electrons along each stream tube. 
This distinction is basic since there exists no a priori theoretical 
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relationship which can suggest how such a "cross tube" sample should appear 
in the fits previously presented. 

Our working hypothesis is to test the idea of SO-II that thermal 
electrons, along the arc length A of the stream tube labeled by r, obey tho 
hypothesis log T(s,A)) = a log n (s,A) + log Z(s) , where a is a constant 
independent of stream tube, but S(s) may change in some way across tubes as 
boundary conditions prescribe. From this relation it follows that 

d d d 

log T(s(t) ,A(t)) = a log n(s(t) ,A(t)) + log S(s(t)) (7) 

dt dt dt 

where the time derivatives are those determined by the time series 
encountered by the detector as provided by the combined motion of the 
spacecraft and the bulk flow of the medium between I .38 and 4. 76 AU. The 
last expression of ( 7 ) vanishes only when the above motion carries the 
detector along the stream tube^ in the flow, and certainly does not vanish 
in any average sense over the protracted time series considered here . 
Rearrangement of equation (7) and noting the properties of total 
derivatives we obtain 

d log T(s(t),A(t)) _ ^ d log r(s(t)) 
d log n(s(t),A(t)) d log n(s(t),A(t)) 

where, as before, the last term is in general non-vanishing. 

The fitting procedure described in the previous section "characterizes" 
the average value of the left-hand side of ( 8 ) . Let us denote this 
procedure by the angular brackets < >; we thus identify 


fit 


^ < 


d log T(s(t) ,A(t)) 
d log n(s(t) ,A(t)) 


> = . . < 


d log 2 (s(t)) 
d log n(s(t) ,A(t)) 


= a + 3 (a) . ( 9 ) 


The non-vanishing size of the last term in (9), which we have defined as 3, 
complicates the simple identification of with the physically 

interesting quantity a of our hypothesis. The magnitude of 3 must be 
assessed before the size of a or its precision may be determined. 

The quantity 3 is not only a function of n and T but also the 

i/' 

physical variable. We initially assume that a ^ to gauge the size of 

3 (a) 3(a^j^^); however, if 3 makes a sizeable correction to this may 
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not be an accurate approximation. As shown in Table 1 the value of 
<P(c»|.it)> (as determined by minimizing perpendicular dispersion of linear 
hypothesis relating log 2 and log n) for all the data used in Figure 1 

is 30/t of We thus conclude that of Figure 1 is not exactly the 

best v?lue to etermine the polytrope index. 

Among the ^.hysical correlations involved in determining the average 
size of 3(a) in our sample are at least the following: 



d log S(a,s(t)) 
d log n(s(t) ,«,(t) 



3 log 2 
3 log u 

3 log 2 
3 log a 


3 log u 


3 log n ' 

(10) 

3 log a 

3 log u 

IL 


/ +. • • T 

3 log U 

3 log n ' 


which says that 3 may be comprised of non-zero average dependences of 2 on 
flow velocity u, density n, or with polytrope index via its bulk speed 
dependence. Since u and n are strongly correlated (Neugebauer and Snyder 
(1966); Burlaga and Ogilvie (1970)), it is unlikely in any given data 
collection that 3 log u/3 log n averages to zero. The next two entries in 
Table 1 shovr the values of a^^^, I, and 3(a^^^) for two disjoint speed 
subsets of the composite data: u ^ 400 km/s. The a fit values of these 

subgroups differ from that of the composite group by 'i‘\% whereas the 3 
values have been reduced more than 200$. The shown in Table 1 have 

also been significantly reduced. However, the average stream tube 
constants 2 = 10^ have also changed significantly and there is a clear 
distinction between 2(u > 400) , and Z(u < 400) with 2 larger in the high 
speed flow. (The errors quoted in Table 1 are 3cr errors.) This is also 
consistent with (2 j^^^) > (^yQygggj,) since the Mariner 10 epoch had a 
higher average speed than the Voyager interval (cf . SO-II and Sittler jgt 
al . . 1979 ). The a fit values of the two speed subgroups are identical even 
at the la confidence level. We therefore conclude that 


that 


3 log 2 

is > 0 

3 log u 

3 log a 
0 


( 11 ) 


9 log U 


(12) 
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and that the complications implicit in (10) are somewhat alleviated by 
narrowino the speed dispersion of the sample. 

The results of equation (11) and (12) and the density-bulK speed anti- 
correlation noted above imply that should be expected to be negative 

provided there are no unexplored functional dependences in equation (10) 

(+ ) . The negative fib values of 3 are therefore compatible with the 

initial fits which determine and I for the two speed subgroups. It is 

worth emphasizing that equation (11) implies that the entropy-like quantity 
Z is positively correlated with the local bulk velocity. 

Fortunately the data set u < MOO kra/s has a very small systematic 
residual, and in this case a '^0.175. We use 

this to refine our estimates of 6 for all data subgroups since 3 is a 
function of the physical ot. This has been done L' the ninth column of 
Table 2. The revised values of a determined in this manner are given in 
the eleventh column. This last entry has been omitted for the first row 
since the composite data set has the largest systematic 3 and the worst 
normalized x* of all the entries in this table. The second and third rows 
of this table yield a weighted average a of a = 0.185 ± 0.015 (Icr) ± 0.045 
(3a) where we have assigned the one sigma uncertainty to correspond to the 
(±) offsets of a required to attain the mean values of data subsets (2) or 
(3). 


The fourth row in Table 1 illustrates that the above errors are 
reasonable by renf.ating our procedure for the most frequently observed bulk 
speed interval (350-400 km/sec) with an estimated physical a 2 ^Q_j^QQ = 

0.205 ± 0.034 (3a), which is almost completely contained in the (3a) 
confidence interval of a given above. 

We estimate that the most probable value is a = 0.175 - 0.030 (So’) 
based on small systematic correction required in data grouping (2) of Table 
1; we have determined a (3o) confidence interval to not exclude any 
defensible estimate of a. In view of the systematic corrections required 
in groups (3) and (4) we do not believe the present data support any 
significant variation of a with bulk speed . 


SUMMARY 

In summary we have come to the following conclusions: 

( 1 ) The experimental results presented here unambiguously support the 
widespread impression in the literature (e.g., Burlaga ^ al . (1971); 
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Hundhausen and Monbgomery (1971); Hundhausen (1973), that the high electron 
mobility in the solar wind allows electrons bo behave more like an 
isothermal than adiabatic gas; 

(2) Our results give further experimental support bo the arguments in 
SO-II that 7 kT^(r) is nearly the local interplanetary potential eoCr) , 
since a pvjlytrope law with a s 1/6 was a necessary and sufficient condition 
for this result. In Figure 3 we graphically exhibit the confidence windows 
for the various suggestions of a made to date. The 99. 9^ confidence 
intervals of the Voyager fib data of the current work does have a 
sizeable overlap with the corresponding 99 *9^ confidence intervale of “ 
suggested by Soudder and Olbert (1979b). 

(3) Indirectly our results also imply that the proper frame electric 
field beyond at least 0.45 AU is dominated by the polarization berm. 

(4) As also shown in Figure (3) the Voyager-Mariner 10 empirical 
results are incompatible at greater than the 99 - 9 % level of confide v.'ith 
a 99. 9? confidence interval of a = .45 recently reported by Feldman ec al . 
( 1978 ). We attribute this difference bo the possible biased sample of 
stream tube constants (d log S/d log n) encountered in that limited study 
and/or the smaller dynamic range of the variables available to the 1 AU 
bound observer. (In the IMP study, the dynamic range of the density was a 
factor of five (5) and only high speed flows were considered. Note that 
the higher speed Voyager subset of Table 1 is also completely incompatible 
with the IMP results) . 

(5) Our composite data set indicates that on a statistical basis the 

one e-folding of the stream tube constant S is found to occur at 50^ of 
the typical value with larger stream tube constantii being more typical of 
higher speed flows. This sense and magnitude of stream tube constant 
variation was that suggested in SO-II to reconcile Y = 7/6 of theory with 
Y “ 1.45 reported by Feldman (1978). 

(6) The suggested value of y = 7/6 made in SO-II is most sensitive to 
an accurate description of the lowest two even moments . n , P , of the 
electron velocity distribution. These moments are dominated respectively 
at approximately the 95 and 75^ levels by the contributions from the 
thermal regime electrons. It is inconsistent to infer from our results 
that the heat flow which is an odd moment over all the electrons is an 
enthalpy flux. The heat flux is controlled microscopically by the 
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asynunetrical distribution in velocity space of the trans- and extrathermal 
electrons (SO-I) which comprise a small fraction of the density and 
pressure. Therefor(s, this empirical corroboration of the polytrope law for 
thermals sheds no direct light on the control of the peculiar, non-dominant 
fraction of the electron phase space which actually implements the heat 
flow of the medium. Rather, these data highlight the control the thermal 
electrons have on the polarization potential which in turn is an important 
force in the subsequent balance of forces which determines the heat flow. 
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FIGURE CAPTIONS 

FIGURE 1 Main figure depicts Voyager data (+) binned as described in 

the text and the optimal power law hypothesis they determine. 
The dots depict all Mariner 10 data as binned in relation to 
the Voyager relation . Error bars are mathematically 
perpendicular to the optimal hypothesis line since the 
averages and dispersions implied are those for the data found 
in bins along the best fit line. 

(Inset 1A) depicts the x* surface for the non-linear 
optimization which selected the coefficients of the line 
plotted in the main figure. This inset emphasizes that the 
local minimum sele.cted is not topologically connected by a 
descending trough to other regimes in the parameter space 
even though the non-linear procedure does not guarantee that 
this need be true. 

(Inset IB) illustrates in the square box the actual 
probability distribution of individual hourly averaged 
Voyager points as a function of the normal distance (scaled 
by the variance within each bin, D » from the best fit 
line selected in the main figure. The right hand vertical 
strip denotes the composite probability distribution of 
observations taken as a set regardless of bin location along 
the best fit line. The upper horizontal code exhibits the 16 
level intensities used in the probability plots. The lower 
horizontal code illustrates the percentage distribution of 
the analyzed data in each bin. Three bins of inset IB are 
used for each point in the main figure. Refer to the text 
for further details and interpretation. 

FIGURE 2 The histogram illustrates the percent relative occurrence of 

stream tube constants Z as a function of departures from the 
average value. The dashed curve depicts a gaussian 
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distribution of equal variance; the observed distribution has 
an e-folding value at 50$ in excess of the typical value 
and is nearly, but not precisely, skew free. This feature of 
the data time series complicates statements of precision as 
discussed extensively in the text. 

FIGURE •R A comparison of confidence fit intervals as suggested by 1) 

the present in situ analysis between 0.45 and 4.76 AU, 2) the 
theoretical suggestion of Scudder and Olbert 1979b (SO-II ) , 
and 3) the recent in situ analysis at 1 AU by Feldman et al . 
(1978) . With the aid of this figure a comparison can be made 
of the agreement or not that is manifest between alternate 
lines of reasoning. Above the 60$ confidence level the 
present analysis and the regime suggested in SO-II are 
indistinguishable; at all but the highest levels of 
confidence the present results are incompatible with those 
suggested by Feldman ^ (1978). 
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